AD717759 


»BBs-tik7  i- o  0  8  ■> 

FINAL  TECHNICAL  REPORT 


A  FOURIER  ANALYSIS  OF  THE  NEAR  FIELD 
VELOCITIES  RESULTING  FROM  A  SIMULATED 
BILATERAL  SHEAR  FRACTURE 

BY 

Merle  E.  Hanson,  Allan  R.  Sanford  and  Ronald  J.  Shaffer 

January  15,  1971 

ARPA  Order  No.:  1031  Amendment  No.  4 

Program  Code  No. :  0F10 

Name  of  Contractor:  New  Mexico  Institute  of  Mining  and  Technology 
Effective  Date  of  Contract:  18  November  1969 

Contract  Expiration  Date:  17  November  1970 

Amount  of  Contract:  $35,000.00 

Contract  No.:  F44620- 70-C-0055 

Principal  Investigator:  Merle  E.  Hanson 

Phone:  (505)  835-5312 

Project  Scientist:  William  J.  Best 

Phone:  (202)  694-5457 

Title  of  Work:  Artificial  Stimulation  of  Earthquakes 


SPONSORED  BY 

ADVANCED  RESEARCH  PROJECTS  AGENCY 


ARPA  ORDER  NO.  1031 

lUTIONALT&SlNICAL  ] 
INFORMATION  SERVICE 

MUM  V*.  331SI  . 

y;  hne  been  e-eo-rovea  for  pucliO 

release  and  sale  ;  its  distribution  is  uninits 


BEST 

AVAILABLE  COPY 


SUMMARY 
of  the 

.  FINAL  TECHNICAL  REPORT 
on  the 
CONTRACT 

ARTIFICIAL  STIMULATION  OF  EARTHQUAKES 
Contract  No.  F44620- 70-C-0055 


A  FOURIER  ANALYSIS  OF  THE  NEAR  FIELD  VELOCITIES  RESULTING 
FROM  A  SIMULATED  BILATERAL  SHEAR  FRACTURE 


The  purpose  of  this  project  is  to  analyze  the  dynamic 
stress  fields  in  geologic  media  particularly  those  created  by 
fracture  or  other  dynamic  source  functions.  Lack  of  realistic 
analytic  descriptions  has  led  to  the  study  of  dynamic  fracture 
in  strained  elastic  media.  An  added  goal  is  to  aid  the 
observational  seismologist  in  the  interpretation  of  transient 
data  from  real  seismic  events  for  the  description  of  source 
type  and  source  parameters. 

A  two  dimensional  simulation  of  a  dynamic  shear  fracture 
was  performed  and  reported  earlier  in  the  contract  year.  The 
velocities  at  several  points  in  the  near  field  have  been  Fourier 
analyzed  to  determine  the  amplitude  characteristics  resulting 
from  the  bilateral  shear  fracture. 

The  holes  in  the  amplitude  spectrums  of  the  velocity  data 
in  the  near  field  are  compared  with  an  analytically  calculated 
source  time.  The  periods  resulting  from  the  holes  in  the 
amplitude  spectrums  correlate  with  the  source  time  calculated 
from  analytical  considerations.  There  is,  however,  the  added 
complications  of  the  near  field  which  results  in  both  of  the 
body  wave  speeds  affecting  some  of  the  points.  The  radiation 
from  either  side  of  the  bilateral  fracture  can  affect  the 
points  being  observed.  The  calculated  time  trace  was  not  of 
sufficient  length  to  indicate  whether  the  complete  spectrum 
was  determined.  The  study  indicates  that  the  spectrum 
analysis  of  the  near  field  velocities  can  show  at  least  part 
of  the  characteristics  of  the  finite  source. 
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ABSTRACT 

The  data  from  a  numerically  simulated  dynamic  bilateral 
shear  fracture  was  analyzed  using  Fourier  spectral  techniques. 
The  holes  in  the  amplitude  spectrums  of  the  velocity  data  in  the 
near  field  are  compared  with  an  analytically  calculated  source 
time.  The  periods  resulting  from  the  holes  in  the  amplitude 
spectrum  correlate  with  the  source  time  calculated  from 
analytical  considerations.  The  analysis  of  the  near  field 
velocities  at  an  observation  point  can  show  at  least  part  of  the 
characteristics  of  the  finite  source. 


INTRODUCTION 


Several  authors,  e.g.,  Davies  and  Smith  (1968),  have  shown 
that  amplitude  properties  of  body  waves  can  be  used  to  determine 
source  parameters  such  as  fault  area  and  fracture  velocity. 

Most  of  these  studies  have  been  applied  to  data  retrieved  at 
distances  from  the  source  that  are  large  compared  to  the 
dimensions  of  the  source.  In  addition,  they  have  been 
primarily  concerned  with  frequencies  lying  within  a  narrow 
band  width. 

In  our  research,  amplitude  and  phase  spectra  have  been 
obtained  f'  r  the  near  field  velocities  associated  with  a 
simulated  bilateral  shear  fracture.  A  discussion  of  some  of 
the  characteristics  of  the  dynamic  shear  fracture  are  given 
by  Hanson,  Sanford,  and  Shaffer  (1971).  The  apparent 
duration  of  fracture  at  a  point  away  from  the  fracture  (hereafter 
called  source  time)  had  to  be  defined  for  radial  distances  from 
the  fracture  comparable  to  the  fracture  length.  Two- space 

dimensions  were  considered  with  the  particle  velocities 
determined  perpendicular  (v)  and  parallel  (u)  to  the  plane 
of  the  fracture. 

TECHNIQUE 

The  simulated  fracture  described  by  Hanson,  Sanford, 
and  Shaffer  (1971)  formed  in  a  block  of  brittle  elastic 
material  infinite  in  the  z  direction  (plane  strain).  The 
block  was  initially  subjected  to  shearing  and  normal  stresses. 
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A  fracture  nucleus  formed  at  the  center  of  the  block  (x  =  0, 
y  =  0)  at  t  -  0  and  propagated  bilaterally  along  the  x-axis. 

The  normal  stresses  were  adjusted  so  that  the  fracture  axis 
was  an  axis  of  maximum  shear  stress.  The  fracture  was 
made  to  travel  at  constant  velocity,  0.4  of  the  P-wave 
velocity,  for  a  given  length  and  stop.  The  final  half¬ 
crack  length  chosen  was  7  centimeters.  The  Lame  constants  and 
the  density  chosen  resulted  in  a  P-wave  velocity  of  0.33333 
cm/microsecond  and  an  S-wave  speed  of  0.19245  cm/microsecond. 
Fracture  acceleration  and  deceleration  were  not  included  in 
the  model  study.  The  calculation  was  terminated  before 
reflections  from  mesh  boundaries  would  affect  the  velocity 
histories  of  the  points  which  were  spectrum  analyzed. 

A  two-dimensional  dynamic  velocity  data  was  generated 
in  the  two-directions,  parallel  (u)  and  perpendicular  (v) 
to  the  fracture,  for  all  the  mesh  points  in  the  grid  by  the  model. 
This  data  was  stored  in  digital  form  with  uniform  time  spacing. 
Examination  of  the  figures  2-11  shows  that  the  velocity  time 
series  were  213.75  microseconds  long  with  171  discrete  data 
samples  per  trace.  A  digital  form  of  the  Fourier  trahsform  _ 

A 00  -  I  a„  cos  (%-)  At  -  il  a  sin  C-jr31)11  CD 

n-=0  n  m  n-0  n  m 

k  ■  0,  1,  2,  . .  m 

was  used,  where  N  represents  the  total  number  of  data  samples 
per  trace  and  m  was  chosen  so  that  the  maximum  frequency  on  the 
figures  is  about  1/2  the  Nyquist  frequency.  The  frequency  is 
0.00117  k  cycles/microsecond  where  k  is  defined  above. 
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Theoretical  studies  (Ben-Menahem,  1962)  indicate  that 
the  velocity  spectrum  should  include  a  factor  of  the  form 
(sin  sin  TrfT)/iTfT  approximately.  This  factor  is  an 
indication  of  the  finiteness  of  the  size  of  the  source  and 
will  produce  holes  in  the  amplitude  spectrum  when  fT  =  n, 
or  for  the  first  harmonic  when  the  period  t  =  T.  For  the  near 
field  the  form  taken  here  for  an  indicated  .source  time  T 
is  taken  as 
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where  is  the  source  time  (as  measured  from  the  right 

hand  tip  of  the  bilateral  fracture)  for  the  right  half  of  the 
fracture  and  T^  is  the  source  time  for  the  left  half  of 
the  fracture.  Tf  is  the  actual  source  time  from  start  of 
fracture  propagation  to  completion  oi  fracture,  RQ  is  the 
radial  distance  from  the  center  of  the  fracture  to  the  point, 
0  is  the  angle  measured  from  the  fracture  axis  to  the  radius 
Rq,  and  b  is  the  half  crack  length.  and  are  body 

wave  speeds  with  the  superscripts  1  and  2  indicating  a  path 
from  the  center  and  ends  of  the  crack,  respectively. 


(2) 
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Normally,  separate  amplitude  spectra  are  obtained  for 
P  and  S  waves.  This  means  that  only  a  single  velocity  need 
be  considered  in  the  interpretation  of  the  spectra.  In  this 
report,  the  spectra  include  both  P  and  S  waves.  Therefore, 
the  interpretation  of  the  holes  becomes  very  much  more 
complicated  because  of  the  interaction  of  P  and  S  waves. 

For  our  spectra,  there  are  four  rather  than  two  possible 
source  times  to  consider  for  each  spectra. 

The  points  analyzed  by  spectral  techniques  were  chosen 
to  be  only  in  the  upper  right,  hand  quadrant.  The  point  of 
fracture  nucleation  is  the  origin  of  the  coordinate  system 
and  the  fracture  lies  along  the  x  axis. 

DISCUSSION  OF  RESULTS 

Shown  in  this  paper  are  the  two  velocities  and  their 
amplitude  spectrums  for  5  points  spaced  about  the  crack  in 
the  upper  right  hand  quadrant.  The  observation  points  were 
chosen  so  that  the  angle  0  measured  from  the  positive  x  axis 
varied  from  19  to  90  degrees.  Figure  1  shows  the  fracture  and 
the  observation  points  approximately  to  scale  for  the  two 
dimensional  simulation. 

Figures  2-11  show  the  velocity  in  the  u  and  v  directions 
as  a  function  of  time  and  the  associated  amplitude  spectrum 
plots  for  the  five  observation  points  near  the  bilateral 
shear  crack.  Table  1  gives  the  angle  between  the  x-axis  and 
the  line  from  the  center  of  the  fracture  to  the  point,  the 
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DISCUSSION  OF  RESULTS,  Continued 

associated  radius  from  the  fracture  center  and  the  source 
time  from  equation  2  best  fitting  the  first  observed  hole 
in  the  spectra.  The  last  column  in  Table  i  indicates 
whether  the  hole  arises  from  movement  on  the  left  side  or 
the  right  side  of  the  fracture. 

Figures  2  and  3  show  the  velocity  and  amplitudes  for 

Rq  =  12.5  cm  and  0  =  90°.  All  the  amplitude  spectra 

indicate  significant  movement  at  the  lower  frequencies. 

This  may  be  partly  due  to  the  truncation  of  the  time  series 

used.  However,  the  time  traces  of  the  particle  velocities 

indicate  plenty  of  low  frequency  motion.  The  parallel 

velocity  trace  shows  a  significant  hole  at  K  *  15 

corresponding  to  a  period  of  57  microseconds  (see  Table  1) . 

r  ii 

Equation  2  (using  the  shear  wave  speed  for  *  and  the 
compressional  wave  speed  for  results  in  a  time  T  of  58 

microseconds.  The  direction  of  the  shear  field  is  toward 
the  right  above  the  fracture.  It  is  possible  that  there  could 
also  be  a  hole  at  lower  frequency  resulting  from  the  exclusive 
use  of  shear  wave  speeds  in  the  equation.  In  this  case,  however, 
the  time  series  was  not  of  sufficient  length  to  show  this  hole, 
v  at  this  point  does  not  exhibit  any  strong  holes  in  its 
amplitude  spectrum.  There  is  a  weak  hole  at  K  *  30  corresponding 
to  a  period  of  29  microseconds.  This  corresponds  to  shear 
wave  body  speed  from  the  fracture  centf-'*  and  compressional 
body  wave  speed  from  both  of  the  stopped  tips. 
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An  example  of  where  a  hole  is  not  well  defined  is  the 
spectra  of  u  at  the  point  3  (Rq  =  14.3  and  0  =  61°)  Figures 
6  and  7.  From  Table  1,  it  can  be  observed  that  the  correlation 
of  the  hole  at  K  =  35.3,  a  period  of  24  microseconds,  is  not 
good  with  the  closest  value  of  T  calculated  from  Equation  2. 

In  addition  the  hole  in  the  spectrum  is  weak.  For  this 
point,  v  exhibits  a  fairly  well  defined  hole  in  its  amplitude 
spectrum  at  K  =  14,  at  a  period  of  61  microseconds.  This 
corresponds  to  the  time  T  (65  microseconds)  obtained  from 
Equation  2  when  the  compress ional  body  wave  speeds  are  used 
for  both  C™  and  C^'.  The  hole  in  the  spectrum  in  this 
case  results  from  the  radiation  from  the  left  half  of  the 
fracture . 

The  point  2  (RQ  ■  22.7  and  0  »  33°)  shows  two  holes 
in  the  spectrum  of  1  il  at  K  =  12  and  25.2,  which  corresponds 
to  periods  of  71  and  34  microseconds,  respectively  (see 
Figures  8  and  9).  Using  Equation  2,  it  can  be  seen  from 
Table  1  that  the  first  hole  results  from  radiation  from  the 
left  side  of  the  fracture  and  the  second  hole  results  from 
radiation  from  the  right  side  of  the  fracture.  In  both 
cases,  the  compressional  wave  speed  was  used  for  and  • 

v  exhibits  a  similar  pattern  for  similar  reasons. 

On  the  other  hand,  it  can  be  seen  fro  .  Table  1  that  the 
first  hole  is  the  result  of  shear  wave  radiation  from  the 
left  half  of  the  fracture. 
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DISCUSSION  OF  RESULTS,  Continued 

Similar  discussions  can  be  made  for  the  other  points. 
Results  are  given  in  Table  1.  In  general,  u  directly  above 
the  right  tip  does  not  exhibit  any  well  defined  holes.  This 
probably  occurs  because  of  the  complicated  overlapping  of 
the  radiation  patterns . 

Complete  interpretations  of  the  spectrum  have  not  been 
made,  primarily  because  a  longer  time  series  should  be  used 
to  determine  the  amplitude  spectrums.  This  procedure  was 
not  possible  in  this  case  as  a  larger  computer  is  required 
to  simulate  the  dynamic  phenomena.  However,  the  amplitude 
spectrum  of  the  near  field  velocities  presented  here  will 
show  at  least  part  of  the  character  of  the  finite  source. 

CONCLUSIONS 

The  analysis  of  the  amplitude  spectrum  shows  that  the 
spectrum  of  the  near  field  velocities  at  an  observation 
point  can  show  at  least  part  of  the  characteristics  of  the 
finite  source.  The  time  series  generated  by  the  simulated 
two-dimensional  dynamic  shear  fracture  was  not  of 
sufficient  length  to  show  conclusively  that  there  are  not 
holes  in  the  spectrum  at  longer  periods.  As  can  be  seen  from 
the  graphs  of  particle  velocity,  the  duration  of  the  time 
series  is  approximately  the  same  length  as  the  fundamental 
wavelength.  A  longer  time  series  would  reduce  some  of  the 
erroneous  amplitude  information  arising  from  tnncation 
of  the  time  series. 
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Application  of  the  spectrum  techniques  to  near  field 
data  with  known  simulated  fracture  phenomena  may  provide 
a  method  of  determining  source  parameters  for  real  seismic 
events.  However,  longer  time  series  (a  greater  number  of 
discrete  data  points)  are  required.  This  imposes  a  restriction 
that  the  simulation  will  require  a  larger  grid,  hence  larger 
computers  and  longer  computer  times. 
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FIGURE  CAPTIONS 

FIGURE  1- -Sketch  approximately  to  scale  of  the  geometry  of 

the  fracture  and  the  observation  points  whose  velocities 
were  Fourier  analyzed.  The  initial  shear  field  and  the 
angle  0  used  in  Equation  2  are  displayed. 

FIGURE  2- -The  parallel  velocity  (u)  and  the  associated 

amplitude  spectrum  for  point  5  on  Figure  1,  the  time 
on  the  velocity  plot  is  in  microseconds. 

FIGURE  3- -The  perpendicular  velocity  (v)  and  the  associated 
amplitude  spectrum  for  point  5  on  Figure  1,  0  the 

time  on  the  velocity  plot  is  in  microseconds. 

% 

FIGURE  4- -The  parallel  velocity  (u)  and  the  associated 

amplitude  spectrum  for  point  4  on  Figure  1,  the  time 
on  the  velocity  plot  is  in  microseconds. 

FIGURE  5- -The  perpendicular  velocity  (v)  and  the  associated 
amplitude  spectrum  for  point  4  on  Figure  1,  the  time 
on  the  velocity  plot  is  in  microseconds. 

FIGURE  6--The  parallel  velocity  (u)  and  the  associated 

amplitude  spectrum  for  point  3  on  Figure  1,  the  time 
on  the  velocity  plot  is  in  microseconds. 

FIGURE  7- -The  perpendicular  velocity  (v)  and  the  associated 
amplitude  spectrum  for  point  3  on  Figure  1,  the  time 
on  the  velocity  plot  is  in  microseconds. 

FIGURE  8--The  parallel  velocity  (u)  and  the  associated 

amplitude  spectrum  for  point  2  on  Figure  1,  the  time 
on  the  velocity  plot  is  in  microseconds. 


FIGURE  CAPTIONS,  Continued 


FIGURE  9- -The  perpendicular  velocity  (v)  and  the  associated 
amplitude  spectrum  for  point  2  on  Figure  1,  the  time 
on  the  velocity  plot  is  in  microseconds. 

FIGURE  10- -The  parallel  velocity  (u)  and  the  associated 

amplitude  spectrum  for  point  1  on  Figure  1,  the  time 
on  the  velocity  plot  is  in  microseconds.  The  high 
frequency  oscillation  seen  on  the  velocity  is  the 
calculational  grid  noise  from  the  numerical  model. 

FIGURE  11- -The  perpendicular  velocity  (v)  and  the 

associated  amplitude  spectrum  for  point  1  on  Figure  1 
the  time  on  the  velocity  plot  is  in  microseconds. 
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data  in  the  near  field  are  compared  with  an  analytically  calculated  source  time.  The 
periods  resulting  from  the  holes  in  the  amplitude  spectrum  correlate  with  the  source 
time  calculated  from  analytical  considerations.  The  analysis  of  the  hear  field  veloc- 
ties'  at  an  observation  point  can  show  at  least  part  Of  the  characteristics  of  the 
finite  source./ \ 
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